Thermodynamics of the one-dimensional half-filled 
Hubbard model in the spin-disordered regime 

S Ejima^ F H L Essler^ and F Gebhard^ 

^ Fachbereich Physik, Philipps-Univcrsitat Marburg, Germany 

^ Rudolf Peierls Centre for Theoretical Physics, Oxford University, 1 Keble Road, 

Oxford 0X1 3NP, UK 

E-mail: f lor ian . gebhardSphysik . uni-marburg . de 

Abstract. We analyze the Thermodynamic Bethe Ansatz equations of the one- 
dimensional half-filled Hubbard model in the "spin-disordered regime" , which is 
characterized by the temperature being much larger than the magnetic energy scale 
but small compared to the Mott-Hubbard gap. In this regime the thermodynamics of 
the Hubbard model can be thought of in terms of gapped charged excitations with an 
effective dispersion and spin degrees of freedom that only contribute entropically. In 
particular, the internal energy and the effective dispersion become essentially indepen- 
dent of temperature. An interpretation of this regime in terms of a putative interacting- 
electron system at zero temperature leads to a metal-insulator transition at a finite 
interaction strength above which the gap opens linearly. We relate these observations 
to studies of the Mott-Hubbard transition in the limit of infinite dimensions. 
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1. Introduction 

The Mott-Hubbard metal-insulator transition in the one-band Hubbard model at half 
band- filling continues to pose an intriguing problem 51121! even in the limit of infinite 
dimensions jSlIll- Starting from the limit of large interactions, the Mott-Hubbard gap 
characterizing the insulating phase diminishes as U is lowered and eventually closes at 
a finite Uc- The determination of the precise value of the latter is an unsolved problem 
in infinite dimensions: The equations from Dynamical Mean-Field Theory need to be 
solved numerically and the best available treatments using the Dynamical Density- 
Matrix Renormalization Group method lead to confiicting results [SHEl- 

Other attempts to determine Uc involve extrapolations of finite-order l/t/-expan- 
sions of the Mott-Hubbard gap [7] and the ground-state energy 0, respectively. These 
expansions can be carried out because the spin background is completely disordered 
and one is left with finding the description of the effective motion of the charge degrees 
of freedom. However, the extrapolated values for Uc obtained by such methods vary 
significantly, depending on the details of the procedure. 

In light of these issues it is desirable to have an example of a Hubbard-type model 
with a disordered spin background which can be solved exactly. One such example is 
the Falicov-Kimball model at half band-filling in the disordered phase . However, its 
ground-state energy turns out to have a trivial l/[/-expansion and hence is of limited 
utility. Curiously, the exactly solvable one-dimensional Hubbard model jTU] at half band- 
filling features a spin-disordered regime at a small but finite temperature in the limit 
of strong Coulomb repulsion, see Refs. [HI [121 CSl • It is then an interesting question to 
analyze the physics of this regime in view of the above issues raised by the studies of 
the Mott-Hubbard insulator in infinite dimensions. 

Our presentation is organized as follows. In section[2]we review the thermodynamics 
of the one- dimensional Hubbard model j^lEinni- In section [HI we analyze the Ther- 
modynamic Bethe Ansatz (TEA) equations in the spin-disordered regime and calculate 
the internal energy and the effective dispersion of the charge degrees of freedom. In 
section [D we interpret these results in terms of a putative interacting-electron system at 
zero temperature. We conclude in section [21 

2. Thermodynamics of the one-dimensional Hubbard model 
2.1. Hubbard model and TEA equations 

The Hubbard model on L lattice sites is described by the Hamiltonian 

H = -tjzil + a;+i,.c,.,J + t/E(%,T - l/2)(%i - 1/2) .(1) 

i=i <T=Ta i=i 

In the following we define u = U/{At) and set t = 1 as our energy unit. Periodic 
boundary conditions apply, and the kinetic energy is diagonal in momentum space with 
the bare dispersion e{k) = —2 cos k. 
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The 'Thermodynamic Bethe Ansatz (TBA) equations' describe the ratios of the 
distributions of hole and particle excitations in thermal equilibrium at temperature T, 
chemical potential yU and magnetic field B eq. (5.60)-(5.65) of [TU]) 
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ln?7„(A) = s*ln(l + ?7„+i)(l + r7„_i)|^ for n > 2 , (3) 

and likewise for ri'^{A) with l/C{k) replaced by C(^)- The convolution operation is 
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2.2. Free- energy density and internal energy 



In terms of the distribution functions ({k) and r7i(A) the free-energy density can be cast 
into the form (see eq. (5.69) of [10]) 

dA;po(A:) ln(l+C(A;))-r / dAao(A) ln(l+r7i(A)) , (6) 

-TT J —OO 

where Cq denotes the ground-state energy of the one- dimensional Hubbard model at half 
band- filling, see eq. (jMj) . Here, the root densities o-q{A) and po{k) are given by 

du Jq{u) exp{—iuA) 
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(71 

+ exp(2M|cu|) ' ^ ^ 

and Jn{x) are nth-order Bessel functions. With the help of © we can transform the 
free-energy density (0) into the form (see eq. (4.21) of ll6j) 

f{T)=C-Ta{T)-Tp{T) , (8) 
dA(To(A)Re 
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With the help of (jH)) the internal energy 

.df{T) 



e{T,u) = f{T)-T- 



dT 



can be expressed as 



3. The spin-disordered limit 

3.1. Definition of the limit 

In the following we consider the TBA equations in the regime 
J<T < A 



(12) 



(13) 



(14) 



for B = and /i = 0, i.e., the half-filled Hubbard chain in zero magnetic field. Here, 
J{U/t — > oo) = At'^/U is the coupling strength for the spin degrees of freedom, and A is 
the gap for charge excitations, A{U/t oo) = U — At. The inequalities ()14|) imply that 



< 1 
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In this limit, the spin degrees of freedom are 'hot', i.e., the charge degrees of freedom 
move in a random spin background. 



e(T) 




> T 



J A 



Figure 1. Internal energy as a function of temperature for the Hubbard model at half 
band-filling and large interactions. 

The qualitative dependence of the internal energy as a function of temperature is 
shown in Fig. Q In the temperature range given by (fT^ the internal energy is essentially 
independent of temperature because the spin degrees of freedom contribute maximally 
to the internal energy whereas there are exponentially few charge excitations. 
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3.1.1. Charge sector We define the dressed energy for cliarge excitations by 

^{k,u,T) = THak)) ■ (16) 

We will see tliat, in tlie spin-disordered regime, tlie dressed energy lias an expansion of 
the form 

The leading term in this expansion is temperature independent, negative and of order 
0{u). Hence 

m = 0{e--/^) « 1 , (18) 
ln(l + l/ak)) = - HCik)) + 0(e-"/^) . (19) 
Equation leads to a significant simplification of the TEA equations. 

3.1.2. k-A Strings Neglecting terms of order 0{e~'^^^), the TEA equations for k-A 
strings take the form 

lnr^;(A)= s*Hl + r^'^^,)il + r^'„_,)l , (20) 

where we define rj^^A) = 0. The solution to pn|l under the boundary conditions (jSj) 
with yu = is [H] 

^;(A) = (^+1)^-1- (21) 

The corrections to rj' are of order exp(^exp(— -u/T)^. 

3.1.3. A Strings Neglecting terms of order 0{e~^^^), the TEA equations for A strings 
take the form 

lnr/„(A) = s * ln(l + r/„+i)(l + Vn~i)\^ , (22) 
where we define 



VoiA) = exp (^--RevO^j - 1 . (23) 

The solution to under the boundary conditions (0) with B = can be obtained by 
iterative linearization, as shown in Ref. ^Hl- The starting point is the solution r]^°^ to 
the equations with r/o = 0, i.e., 

r^i°)(A) = (n + 1)2-1. (24) 

One then linearizes ()22|) around the solution ()24|) by writing //„ = 77^°^ + r]l^^ (for n > 1) 
and keeping only the terms linear in ri^\ This gives the following set of linear integral 
equations 
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The set of hnear integral equations ()25|) can be solved by Fourier transform ^Hl- We 
find 
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We see that we have 
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as required. 



(27) 
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3.2. Dressed energy of the charge degrees of freedom 
Combining the results for rii{A) and ri[{A) we obtain 

'i+v[{Ay 



In 



Substituting this back into (0) we obtain the first term of the expansion (fTTj) 



K(o)[k,u) = -2cosk-2u- - I —Ji{uj)e 
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(29) 



(30) 



2 J-oo UJ 

Crucially, the contributions of r]i and r][ to ({k) do not feed back into the integral 
equations for ?7„ and r]'„. The reasons for this are that 

(i) the corrections do not change the fact that C is 0{exp{—u/T)) and terms involving 
ln(l + C(^)) can therefore be dropped from the TBA equations; 

(ii) the corrections to depend only on sin(A;) and therefore do not contribute to 



dk cos(/c) s{A — sin k) ln[(^{k)] . 



(31) 



The result flSUjl should be compared to the corresponding expression at zero temperature, 
see eq. (7.10) of fTU], 

Ko{k,u) = -2cosk-2u-2 I — Ji(cj)e^""" \ , . (32) 



f°° du 

-2cosA; — 2u — 2 / — Ji[uj)e 

J-oo CO 

Performing the integral in (jHn|l we finally find {u = U/A) 

1 



1 + exp(2|u;|n) 



fi:(o)(A;, U) = e{k) - 



[e(fc)]2 + [/2 + J[e{kY - t/2]2 + [wy 



(33) 



2^2 

where e{k) = —2 cos A; is the bare dispersion. 

The interesting point is that K^^\k) ^ ^^{k) despite the fact that T -C |fi:o(^)|- 
The reason for this is that spin and charge degrees of freedom are still coupled in the 
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Hubbard model. This is obvious from the TEA equations and also from the known 
scattering matrix of elementary excitations ^Tj, see Chap. 7.4 in ^U]. We note that 
the modification of the dressed energy of the charge degrees of freedom by the spin 
sector is a general characteristic of the spin-disordered regime. In the asymptotic low- 
energy regime, spin and charge degrees of freedom decouple quite generically for one- 
dimensional interacting-electron systems. However, on the scale of the temperature in 
the spin- disordered regime such a decoupling no longer holds. This in turn is expected 
to lead to modifications in non-universal physical properties such as the charge velocity. 
Our result is in agreement with this expectation which we note equally applies to a less 
than half-filled band. 



3.3. Dressed momentum of the charge degrees of freedom 

The total momentum in thermal equilibrium is equal to zero by virtue of translational 
invariance. In order to determine the dressed momentum of our charge excitation, we 
start from the expression for the contribution of a charge excitation with real /c's to the 
total momentum jIH| (see (5.97) and (5.38) of fTU]) 

m = k + j: dAei K^(A) + <(A)] , (34) 

„_i J — oo \ nu I 



where 



d{x) = 2arctan(a;) (35) 
and the root densities for particles and holes obey (see (5.41) of 10 ) 

1 oo /-oo 

ff{k) + p'^ik) = — + COS A; V / dAa„(sin A; - A) [a'P{A) + <(A)] ,(36) 
2n „=i^-oo 
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with 



^"(^) = S7IT ' ^n(^) = ^, ' C(^) = ^- (39) 
an{A) (Tn(A) pP[k) 

As ({k) = (9 (exp(— n/T)) we can drop p'^i^k) from equation (jHBj) and substitute the 
resulting equation into (j37jl and (j38|l . We obtain the following result for the driving 
term 

/TT rTT A Ir 

dks{A - sin k)pP{k) = / — s(A - sin k) = ao(A) . (40) 
-TT J~TT 2TX 
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The temperature-independent contribution to the particle and hole root densities is 
obtained by using (PT|) and ()24|) and then solving the resulting sets of coupled linear 
integral equations. We find 

ar(A) = 0(exp(-VT)) (41) 



and 
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Substituting PHj) back into (jMj) we obtain the temperature-independent contribution 
to the dressed momentum 



p{k) = k -\- - arcsin 



2 sin k 



4^2 + (sin A; + 1)2 + JAu^ + (sin A; - 1) 



(44) 



We note that the physical momenta of holons (li) and anti-holons iK) are obtained from 

p{k) by ph{k) = Tx/2- p{k) = pj^{k) + n. 

3.4- Effective dispersion for strong coupling 

The effective dispersion is given implicitly by equations (j33j) and (j44j) . In order to have 
a consistent expansion for the Hubbard model in the spin-disordered regime, we need 
to truncate these two equations at order 1/u because the next subleading temperature- 
dependent contributions are of the order l/{uT). Hence, the effective dispersion in the 
spin-disordered Hubbard model is given by 



K^^\p,u) = -2u - 2 cosp - ^{3 - 2 cos'^ p) + O (u' 
For the Hubbard model at zero temperature we find instead 



Ko{p,u) 
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Both formulae can be cast into the form 



1 



K(p,u,T) = -2M-2COSP- — (l-47JT))(3-2cos2p) + C(M"2) ,(47) 

4m ^ ^ 



where 



7s(T) = (SiSj+i), 



(4^ 



One- dimensional Hubbard model in the spin- disordered regime 



9 



denotes the nearest-neighbor spin correlation function in the Heisenberg model at 
temperature T. Thermal averages of operators A over spin configurations are defined 
by 



Tr " 


exp (-/^f^Hois) 


A 


Tr 


exp (^-pHueis 


)1 



with the Heisenberg Hamiltonian (see Appendix 2. A of PH] ) 



-f^Heis — X! 



4f 
U 



S,:S 



(49) 



(50) 



In fact, we have 7s(T = 0) = 1/4 — ln(2) from the Bethe-Ansatz solution of the 
Heisenberg model [1^] and %(T = oo) = for uncorrelated spins. Results for all 
temperatures can be found in Refs. I21j. 

The result (j47j) can be obtained within the l/[/-expansion [7j. This approach can 
be used here because the one-dimensional lattice is a Bethe lattice with coordination 
number Z = 2. To leading order in 1/m we must treat the Hamiltonian 

hi = -^PofPifPo , (51) 

where Pq projects onto the subspace of zero double occupancies. At half band-fiUing, 
hi reduces to the Heisenberg model (I^Uj). The internal energy density of the Heisenberg 
model is given by 



es(T) 



1 

u 



7s(T) 



(52) 



and can be calculated analytically in terms of the Thermodynamic Bethe Ansatz P^l^ 
or, equivalently, via the solution of a set of coupled integral equations j2S|- 

For the derivation of the effective dispersion we have to solve (see (43) of [7j) 



U 



9i,o 



= 25-1,2 + 51,0 , (53) 

where ho describes the the hopping of holes in the spin background, Iiq = PqTPq. The 
expectation value on the left-hand-side of (j53|l is readily calculated so that we find 



-8(7s(T)-- 



251,2 + 5i,o 



(54) 



as our first equation. The second equation we obtain from the solution of eq. (44) of [7j, 



U 



^1,2(^^0)^ + ^1,0 {hofc,^ 



= 651,2 + 251,0 . (55) 
The expectation value on the left-hand-side of is readily calculated and we find 



4(7s(r)-- 



25 



1,2 



(56) 
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as our second equation. From this gi^2 = —(1 — 47s(T))/2 and gi^ = 3(1 — 4:%(T)) 
result, and the effective Hamiltonian for the motion of a single hole in a spin background 
becomes 

^„ ^ (-(MVS + 3) , . 

h^'' = ho + {l- 47s(r)) ^ ^ + O {U^') . (57) 

Replacing ho — > — e(p) = 2 cosp and + 2n = —h'^^ as in we find for k{p, u,T) = uj 

Kip, u, T) = -2u-2 cosp - — (1 - 47,(T)) (3-2 cos^ p) + (u-"^) , (58) 

4m ^ ^ 

as used in (jUj). 

From the l/f/-expansion we can determine the density of states of the lower 
Hubbard band. As in the shape- correct ion factor to first order is found to be 

s(e) = 1 - (1 - 47s(r))e/f/ so that we find [«(T) = 1 - 47s(T)] 

4Sb(^) = £ depo(e) (l - a(T)^) sL + U/2 + e + a{T)^ " 



2U 



Po 



U - ^(1 + a(T))f/2 ^ 2a{T)Uu + 6a(T)2 ] /a(T) 



,(59) 



where Po(e) = l/(vr-\/4 — e^) for |e| < 2 is the density of states for non-interacting 
electrons and uj- < uj < ujj^ with uj± = —U/2 ± 2 — a{T)/U. In particular, for the 
single-particle gap we find 

A(i) (u, T) = -2u+ = 4u - 4 + ^ "^^'^^^^ , (60) 

2u 

up to and including the first order in the strong-coupling expansion. 

Finally, we note that the momentum distribution can also be determined along 
these lines. We find 

{n,)s = E(V)s = 1 - ^ (1 - 47s(T)) + O {u-') , (61) 
in agreement with eq. (3.2) of Ref. 



3.5. Internal energy 

We must evaluate (fT^ in the spin-disorder limit. We have 

^ T'— ^ J ^dkpoik)K^'\k,u) , (62) 

e,^T^§; = 0. (63) 

The latter follows from the fact that r][ is independent of temperature. This leads to 
the following expression for the internal energy 

/■°° duj 

e(T,u)= -u- — Jo{iu)Ji{iu)exp{-2uiu) + 0{l/Tu) 

Jo UJ 

= _ (-^\ + 0{l/Tu) = e{u) + 0{l/Tu) . (64) 
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From the internal energy we can derive the average double occupancy, d{u) = 
{l/4){d[e{T,u)+u])/{du), as 

din) = l-^K (-1) + Oil/Tu) . (65) 

Here K(m) and E(m) are the complete elliptic integrals of the first and second kind, 
respectively. The internal energy is to be compared to the ground-state energy of the 
Hubbard model at half band-filling 

r°° duj 1 

eo{u) = -u-4 —J,{uj)Mu;)— — - . (66) 

Jo uj 1 + exp [Zuuj) 

The average double occupancy for the Hubbard model at zero temperature follows from 
the derivative of with respect to the interaction strength. 



4. Zero-temperature interpretation 

4.1. Mott-Hubbard transition 

Now we interpret our results in terms of a putative one- dimensional interacting-electron 
system at zero temperature. In practice, we use the results for the temperature- 
independent contributions to the internal energy and the effective dispersion derived 
for the spin-disorder regime in the Hubbard model, J <C T -C A, for any value of U. 
Our motivation are studies of the Hubbard model in infinite dimensions where the Mott- 
Hubbard insulator is in the spin-disordered phase above the Mott-Hubbard transition. 
In order to model such a situation in a one-dimensional system one would need to take 
the hmit T — after letting J 0. Of course, this is not possible for the Hubbard 
model. 

The one- dimensional Hubbard model at half band-filling describes a Mott-Hubbard 
insulator for all U > 0, i.e., the gap Aq{U)/2 for single-particle charge excitations is 
finite, Ao{U) = -2ko{±h,U). From ^ we have 



Ao(f/) 



-2 + - + A —Muj) 



(67) 



2 Jo UJ ^ ' l + exp{Uuj/2) 

For U —>■ 0, the gap is exponentially small whereas it increases linearly with U for large 
interactions. The transition at Uc = O"*" is readily understood as the consequence of the 
perfect-nesting property in one dimension so that the (marginally) relevant Umklapp 
scattering processes drive the system into the insulating phase for all t/ > 0. 

In our putative interacting-electron system these scattering processes are rendered 
ineffective by the random spin background which, at zero energy cost, provides a 
mechanism to dissipate momentum in scattering processes involving charge degrees of 
freedom. Therefore, we expect that the charge gap will open at a critical interaction 
strength. Indeed, from (jHUj) and we find 



AiU) 



-2 + - + 



/ — Ji(cj)exp(-[/u;/2) 

JO UJ 



-4 + +4 . (68) 
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The gap opens linearly with slope ^/3/2 at 

Uc = = ^ 0.8Q6W , (69) 

where W = 4 is the bandwidth of the Hubbard model. The gap is shown in Fig. |21 




Figure 2. Gap for single-particle excitations as a function of the interaction strength U 
for the interacting-electron system with a disordered spin background at half band- 
filling. 



4-2. Average double occupancy 

We expect that physical quantities such as the average double occupancy display 
unphysical behaviour in the region U < Uc- In fact, d{u) contains a term proportional 
to u ln{u) so that its derivative diverges logarithmically for u 0. This diverging slope 
is seen in Fig. |3] where we compare the average double occupancy of the interacting- 
electron system with a spin-disordered background with the double occupancy of the 
Hubbard model. 



4-3. Strong-coupling expansions 

The internal energy of the Hubbard model at zero temperature and in the spin- 
disordered case can be expanded in powers of 1/U. In both cases, the radius of 
convergence is given by = 4, see (6.83) of [10] for cq^u) and (17.3.12) of ^2^ for 
e{u). Explicitly, 



-2m 
m=l 



[eiU) + U/A]/U = «2mt/ 
=1 

1 



1 3 20 175 1764 19404 , , 



One- dimensional Hubbard model in the spin- disordered regime 



13 
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^^0.15 




0.05 



Figure 3. Average double occupancy do{U) of the Hubbard model at half band-filling 
as compared to d{U) H65|l for the interacting-electron system with a spin-disordered 
background. 



In fl52|) we determined the coefficient 02 = —1 within the l/f/-expansion. From the 
coefficients a2m, one may actually deduce the radius of convergence of the series, as 
done in P, by extrapolating the ratio of the coefficients r(m) = \a2m/(i2m-2\ > 2) 
for m — i> 00. In our case r(m) = 4 (2m — l)(2m — 3)/m^ is a second-order polynomial in 
1/m and we correctly find 

Ul, = ^im (v^) = 4 . (71) 

Note, however, that the radius of convergence of the energy is not related to the critical 
interaction strength of the metal-insulator transition, Uc = 2^/3. Therefore, it may also 
be doubted that this approach ^ is justified for the case of the infinite-dimensional 
Hubbard model. 

In contrast to the internal energy, the series expansion of the gap in the spin- 
disordered Hubbard model converges for U > = 2. The ffist terms of the expansion 
read 

00 r) rj 

A(f/)= E^b^U-"^ = U-A + -^--^±... . (72) 

In (pTn|) we determined the coefficients 6_i = 1, 60 = 4, and bi = 2 within the 1/U- 
expansion. 

Now that Uc = 2v^ is larger than the convergence radius of the series, U^^^ = 2, 
the gap opens linearly and the ffist few orders of the 1 /[/-expansion provide a good 
description of the gap, as shown in Fig. |21 This is particularly true for the critical 
values for the closing of the gap as inferred from the truncated 1 /[/-expansion. Let us 
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denote by A^"^\U) the mth-order truncation of the series, e.g., 
A(°)(f/) =f/-4, 
AW(t/)=f/-4-^, 




(73) 



etc., and let Uj^"^^ be the critical interaction strength at which the mth gap opens, 
= 0. We then find 



The numbers in square brackets give the percentage difference to Uc- It is seen that the 
series converges very fast to the exact value. This observation supports the application 
of this approach to the Hubbard model in infinite dimensions [3IZI- 

5. Conclusions 

We have analyzed the Thermodynamic Bethe Ansatz equations for the half-filled one- 
dimensional Hubbard model in the spin-disordered regime, t'^/U <^ T <^ U. We 
have derived explicit expressions for the leading terms in the internal energy and the 
dressed energy and momentum of the charge degrees of freedom. The resulting effective 
dispersion of holons and anti-holons differs from the corresponding result in the half- 
filled Hubbard model at zero temperature to order J = At^ /U . This effect is due to the 
coupling of charge and spin degrees of freedom and occurs at the expected energy scale. 

We then interpreted the entire temperature-independent part of the effective 
dispersion and the internal energy in terms of a putative interacting-electron system 
at zero temperature. From these results we derived some implications for the analysis 
of large-coupling expansions for the Mott-Hubbard insulator in infinite dimensions. 
Moreover, our analytical result can be used to assess the quality of other approximate 
schemes which describe the effective charge dispersion in a random spin background, 
see, e.g., Ref. pHj . 
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